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Recent progress on the bottomonium states is reported. This talk briefly reviews the observation 
of hb(lP),hb(2P), Z^, Z^ 2 states, transition of hb(nP) — > nb(rnS)^ and new studies on the rjb{^S) 
state. Other n and tv + tv~ transitions of T(nS) are also discussed. 



I. STUDY OF BOTTOMONIUM 

The history of bottomonium starts in mid 1977 where an enhancement was observed at a fi + \i~ mass peak near 
9.5 GeV, from the collision of 400 GeV protons on nuclear targets at Fermilab [1 j. In 1978, with much improved 
resolution, the T particle was observed with a width limited by energy spread of beams at the DORIS e + e~ storage 
ring. Just like the J/ip particle, the T discovery comes from direct production, thanks for the spin-parity J p = 1~ 
and its substantial leptonic partial width. 

However, other particles of the bottomonium family can not be directly produced at an e + e~ collider. So, alter- 
native experimental methods should be used to search for these states. For example, the P-wave states of b and b 
were discovered using the inclusive photon spectrum of Y' decay at Crystal Ball collaboration in 1985. In Fig. [I] left, 
A triplet of peaks between 100 and at 200 MeV is clearly seen, which corresponds to Y' — >• 7X&( 3 -p2,i,o) transitions, 
while the peak at the right is due to reflection from Xb ~^ 7Y. 

Now, a large of set bottomonium states has been predicted and studied, which was nicely arranged according to the 
bb level scheme (PDG). Recent discoveries of new bottomonium states also follow the history, where the inclusive or 
semi- inclusive methods play important role. As two examples not long ago, BaBar collaboration discovered the lowest 
lying bottomonium state rjb using the same inclusive 7 spectrum method [3], and CLEO collaboration discovered 




FIG. 1: Left figure shows the inclusive photon spectrum obtained by the Crystal Ball Collaboration at DORIS II [2 . Right 
figure shows the inclusive missing mass spectrum of 7r + 7r - obtained at the Belle's 121.4 fb _1 T(5S) data [9], where the bottom 
plot is the background subtracted residual distribution, with the fit function overlaid. Different dashed lines indicate the 
position of various resonances. Notice the similarity between left and right figures. 
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T(1D) using the recoil mass of two soft 7s with / + /~ tagging [4 j, which represents the inclusive property of two 7s. In 
the following sections, recent results of bottomonium, especially those obtained from Belle's large 121.4 fb _1 T(5*S) 
data, are discussed. 

II. T(55) AND h b PARTICLES 

In 2008, Belle collaboration reported surprising large partial widths of T(5*S) — » T(l, 2, 3S)tt^7t~ [5 , which are 
two orders of magnitude larger than the widths of T(2,3,46') — >> T(16')7r + 7r _ . While the reason of this discrepancy 
is unclear, an T state different from the conventional bb bound state has been considered [6]. 

In 2011, CLEO-c collaboration suggested that the cross-sections of /i c 7r + 7r~ and J/^7r + 7r~ are of similar magnitude, 
and are all enhanced near 4.26 GeV e + e~ energy [7 , although for h c 7T^7T~ the statistics was limited. Interpreting this 
enhancement as the y(4260) resonance, the partial width of F(4260) — » J/^7r + 7r _ has been shown to be greater than 
508 keV at 90% confidence level [8 , which is much larger than the partial widths of ijj' (102 keV) and (53 keV). 
This behavior is similar to the bottomonium case of T(5S). Thus, we naturally think that the hb7T^7T~ production 
rate is also greatly enhanced in the T(5*S) region, by assuming the similar mechanism as in the charmonium case. 

The hb decay modes are unknown and should be complicated, since it will decay to light quarks. In order to study 
the decay T(55) —> hb7r+7r~, we can still use inclusive method to avoid reconstructing hb directly. The recoil or 
missing mass of 7r + 7r~ " defined as MM(7r+7r-) = y/(P(T(5S)) - P(7r+7r-)) 2 was studied at Belle. Here P is the four 
momentum of relevant particles, with the P(T(5S)) obtained from beam energies. The spectrum of MM(tt^7t~) 
shown in Fig. [I] right is of huge statistics which is around 1 million events per 1 MeV. Nevertheless, it shows the 
similar behavior compared to the Crystal Ball's inclusive photon plot in Fig. [I] left. By fitting it, Belle made the 
first observation of hb(lP) and hb(2P) states [9], whose masses agree with the theoretical expectation from center-of- 
gravity of Xb states. However, the mechanism of T(55) —> hb(nP)7r+7r~ decay is exotic, because the ratio of spin-flip 
to non-spin-flip branching fractions r(T(5£) —> hb(nP)7r+7r~)/r(Y(5S) —> T(26')7r + 7r _ ) are not suppressed. 

III. Z b PARTICLES 

While CLEO-c did not report any resonance structure in the /i c 7r + 7r _ system, Belle has enough statistics to study 
the resonance structure of the three body /i67r + 7r~ in T(5S) decay. This was achieved by looking the missing mass 
of a single 7r + or 7r~, which should effectively be the M{hb^~) or M(hb7r+). Because of symmetry transposing 7r + 
and 7T~, the two missing mass distributions is combined and upper half of the available range is used, which we 
denote as MM(tt) distribution. Then the missing mass of two pions MM(tt^7t~) was fitted to extract hb(lP) and 
hb(2P) signals in bins of MM(tt). The resulting spectra of hb yields as a function MM(tt) will be background- free 
distributions for M{hb{lP)^) and M(hb(2P)7r). The distribution for M{hb{lP)ir) from data exhibits a clear two- 
peak structure without significant non-resonance component (Fig. [2] left). The distribution M(hb(2P)ir) behaves 
similarly with smaller statistics. These two structures are referred as Zb\ and Z52 and parameterized as two P-wave 
Breit-Wigner amplitudes. The fit function with yfs = M(hb7r) is: 

|BW(s, Mi, Ti) + ae^BW(5, M 2 , T 2 ) + be^\ 2 ^L. (1) 

Here the amplitude of two resonances and a non-resonant component are added coherently to form the rate, which 
is then multiplied with a phase-space factor where q or p is the momentum of the pion from T(5S r ) or Zb decay 
in the rest frame of their mother particles. 

The Zb particles can also be studied from the decay T(55) — » T(n6')7r + 7r _ , with n = 1, 2, 3. In this case, the final 
state particles can be fully reconstructed with T(nS) — » /i + \i~ . Dalitz analysis is then performed to the three-body 
final states to extract maximal information. The amplitudes used in the parametrization includes two Breit-Wigners 
for Zbi and Z52, a coupled-channel Breit-Wigner for /o(980) scalar, a Breit-Wigner for / 2 (1270) tensor state, and a 
non-resonant amplitude. Results of Zb parameters for all five channels are consistent, and the average masses and 
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FIG. 2: The M(h b (lP)it) (left) and M[T(25)7r] max (right) distributions in the Z h observation study at Belle [10]. Two peaks 
in both plots show the Z b i and Zt2 signals. 



widths are M = 10607.2 ± 2.0 MeV, T = 18.4 ± 2.4 MeV for Z bl and M = 10652.2 ± 1.5 MeV, T = 11.5 ± 2.2 MeV 
for Z b2 [M- 

If we compare the relative phase of two Z^s in h b (lP)7T^7T~ channel and T(25)7r + 7r _ channel, we found that the 
former is close to 180° and the latter is close to zero degree. From the distribution of M(Y(25)7r) max , as shown in 
Fig. fright, a large destructive interference dip near 10.6 GeV/c 2 suggests that this interference effect exists in the 
whole Dalitz plane, which is only possible if the Z b has a spin-parity J p = 1 + . 



IV. rib PARTICLES 



h b is expected to have large decay branching fractions to 7757 [TT]. The 775 particle previously has been studied [3 , 4 
at BaBar and CLEO in the channel Y(35) —> 7757. Since large h b samples are available at Belle, it is natural to study 
the transition from h b to 7757. 

The method to extract 775 signals at Belle is similar to what is used in extracting h b signals. Here the missing mass 
of three particles MM(7r + 7r _ 7) is calculated. To avoid the correlation between MM(7r + 7r _ 7) and MM(7r + 7r _ ), their 
difference is used, as AMM(7r + 7r"7) = MM(7r + 7r"7) - MM(7r + 7r") + M(h b ), where M(h b ) is fixed at the nominal 
value. The MM(tt^7t~) spectra are then fitted in different AMM(tt^7t~j) bins, to obtain the h b {lP, 2P) yields as a 
function of AMM(7r + 7r _ 7), as shown in Fig. [3] (a,b,c). The fact that h b is produced from the Z b resonances is used 
to suppress the background by the requirement 10.59GeV/c 2 < MM(tt) < 10.67GeV/c 2 . 

The h b (lP : 2P) yields as a function of AMM(7r + 7r~7) in rj b (lS) and rj b (2S) mass regions are fitted to a sum of the 
r] b signal component and a smooth background component, for the transitions h b (lP) — » 775(15)7, h b (2P) — >• 775(15)7, 
and h b (2P) — >• 775(25)7. The signal is described by the convolution of a non-relativistic Breit-Wigner function with 
the resolution. At the time of this report, the fitted masses and widths are m r]b ^ 1 s) = 9402.4 ± 1.5 ± 1.8 MeV/c 2 , 
T mm = 10.8±|;?±|;g MeV, and m Vh{2S) = 9999.0 ± 3.5t 2 j MeV/c 2 [12 . The hyperfine splittings m T(n5) - m m(nS) 
are determined as AM# F (15) = 57.9±2.3 MeV/c 2 , AM HF (2S) = 24.3^;g MeV/c 2 . This agrees with the theoretical 
calculations [T3] , S.Dobbs et al. studied the f] b (2S) in the process T(25) — > 7775(25) in CLEO data, by reconstructing 
26 exclusive modes together for the 775(25) candidate. The hyperfine splitting is reported as AMh F (2S) = 48.7 ± 
2.7 MeV/c 2 [14 . Fig. [3] (d) shows the fit. This is of about 5 sigma discrepancy compared to the Belle result and is 
in strong disagreement with theory. So further experimental clarification is needed. 
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rjb(2S) region is shown in (c). Clear rjb(lS) and rjb(2S) signals are all seen in (a,b,c). The mass difference between T(2S) and 
rjb(2S) candidates in CLEO data [14 is shown in (d), where the unidentified small peak around 50 MeV is interpreted as the 
rjb(2S) signal. 

V. n TRANSITIONS OF T(raS) AND OTHER T(55) DECAYS 

The transition T(nS) — » rjT(mS) is a spin- flip E1M2 transition. From the QCD multipole formalism [15], this 
spin-slip amplitude scales as l/ra^, and this n transition is suppressed compared to the tttt transition. However, the 
experimental values do not support those predictions. Scaling from the known branching fraction ip' tjJ/x/j, the 
branching fraction of T(2S) — >• Y(l S)r) should be around 8 x 10 -4 , but the experimental value is around 2 x 10 -4 [T6] . 
In addition, from Ref. [17 , the relation B(T(4S) T(15>) ~ 2.5 x B(T(4S) T(1S>+tt-) contradicts with the 
expected suppression in n transition. So it is important to study more in these channels. 

Belle performed a study of Y(2S) —> (77, tt°)T(1S) using Belle's 24.7 fb~l T(2S) data. The reconstruction of 
T(1S) — >• e + e~, and a small total reconstructed momentum in the T(5S f ) center-of-mass frame are required. By 

fitting the rj candidate and 77 masses, as shown in Fig. we obtain B(T(2S) -> T(lS)rj) = (3.41 ±0.30 ±0.35) x 10~ 4 
and B(T(2S) T(15)tt ) < 4.3 x 10" 5 (90% CL). 

The T(5S) data was also analyzed similarly to search for transitions to T(l, 25')7r + 7r _ 77 states, with Y(1,2S') 
reconstructed in the channel. If we require the n — >• 7r + 7r _ 7r° selection criterion for the pions and photons, n 

transition of T(55) can be studied. At the same time, if we require the T(25) — > T(15')7r + 7r _ selection, T(55) — > 
T (25)77 can also be studied, since the remaining two 7s can make an 77. This time we simply fit the difference of 
the T(bS) candidate's energy and beam energy in the center-of-mass system to extract the signal. This method 
gave the same result as fitting the missing mass of 77. We found B(T(55) —> T(lS)rj) = (7.3 ± 1.6 ± 0.8) x 10 -4 , 
B(T(55) T(25)t7) = (38 ± 4 ± 5) x 10~ 4 , and B(T(55) T(15)t7 / ) < 1.2 x 10~ 4 (90% CL). Fig. |] (a,b) shows 
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FIG. 4: n (left) and 77 (right) mass distributions in the Y(25) — > nY(lS) and Y(25) — > 7r°Y(15) search at Belle, n candidates 
from 77 and 7T + 7T~7T channels are combined. 




MM(ti), GeV M(yy), GeV MM(iV), GeV 



FIG. 5: Plots of Y(55) — > Y(l, 25)7r + 7r _ 77 study at Belle. Look-back plots for the missing mass of n in the Y(55) — » 
T(l, 25')?7(7r + 7r _ 7r ) study (a) and the 77 mass in the Y(55) — > Y(25)(Y(15)7r + 7r _ )?7 study (b) are shown. Missing mass of 
7r + 7r - after the Xbi,2 — > T(lS)j selection is shown in (c), where the left, middle, right peaks are due to Y(25)7r + 7r _ , the signal 
r(lL>)7r + 7r-, and the X (25)(T (ISV+tO^TT) reflection. 

the look-back plots of the first two modes. 

A peak of T(1D) was seen in the inclusive distribution of missing mass of 7r + 7r - where the states was observed 
(Fig. [T]). A more complete study of the decay Y(55) —> Y(lD)7r + 7r _ is now possible when the exclusive mode 
T(l, 26')7r + 7r _ 77 is being analyzed, since we can aim for the sequential decay of T(1D) — » Xbl —> Y(1S)jj. The 
peak of T(1D) became clearer after this exclusive reconstruction, and turned out to be even sharper after the 
Xbi, 2 T(1*S)7 selection. From the fit shown in Fig.[5](c), we obtain the observation of this channel with £>(Y(55) —> 
T(LD)7r+7r-) x B(T(LD) -+ XbiA lp h) x B(xbiA lp ) r ( is h) = ( 2 -° ± °- 4 ± °- 3 ) x 10 ~ 4 - 

VI. SUMMARY 

Belle's large Y(55) data provided many new results of bottomonium states. The inclusive spectra of 7, 7r°, 7r + 7r~, 
etc. provide the main method used in the discovery of new states from 27 years ago till now. Four new bottomonium 
states /i&(lP), /i&(2P), Z^ 2 were observed at Belle, where the Y(55) behaves like the particle Y"(4260) in the 
charm sector. Observations of the /i&(l,2P) — > 775(15)7 and hb(2P) — >> 775(25)7 were made. The particle 775(25) 
was observed for the first time and the hyperfine splitting AMhf(2S) at Belle agrees with theory. S.Dobbs et al 
obtained a higher value AMhf(2S) from CLEO data and two results are in clear disagreement. For the other 
transitions of the T(55), observations of T(55) — >• Y(1,25)t7 and T(55) — >• T(1D)tt^7t~ channels were made. For 
the 77 transition of other T particles, The new measurement of T(25) — »• T(15)t7 branching fraction is still smaller 
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than theory prediction. Finally, In the talk covered by Andre Chisholm, ATLAS and DO has observed the Xb(3P) 
state from the transition Xb(3P) — >• T(1,2S)7- In future, we expect more Belle results, and new results from the 
Hadron colliders. 
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